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ABSTRACT 

The most general form of the nonrelativistic Grad-Shafranov equation describing anisotropic pres- 
sure effects is formulated within the double adiabatic approximation. It gives a possibility to analyze 
quantitatively how the anisotropic pressure affects the 2D structure of the ideal magnetohydrodynamical 
flows. 

Subject headings: MHD — star: winds, outflows — ISM: jets and outflows 



1. INTRODUCTION 

The Grad-Shafranov (GS) equation for axisymmetric 
stationary ideal magnetohydrodynamical (MHD) flows 
has been considered very intensively within the last two 
decades both for nonrelativistic stellar (solar) wind in- 
cluding outflow from young stellar objects (Okamoto 1975; 
Heinemann & Olbert 1978; Blandford & Payne 1982; Hey- 
vaerts & Norman 1989; Sakurai 1990; Pelletier & Pudritz 
1992; Shu et al 1994), and for relativistic outflows from ra- 
dio pulsars (Ardavan 1979; Bogovalov 1992) and 'central 
engine' in active galactic nuclei (Phinney 1983; Begelman, 
Blandford & Rees 1984; Mobarry & Lovelace 1988; Chi- 
ueh, Li & Begelman 1991; Beskin, Kuznetsova & Rafikov 
1998). As is well known, in this approach the isotropic gas 
pressure has been postulated (see e.g. Heyvaerts 1996 for 
introduction). 

It is necessary to stress that, actually, thermal effects 
were considered rather carefully only within the given 
poloidal magnetic held approximation (Weber & Davis 
1967; Kennel, Fujimura & Okamoto 1983; Paatz & Ca- 
menzind 1996), or within the self-similar approach (Low 
& Tsinganos 1986; Contopoulos & Lovelace 1994; Sauty 
& Tsinganos 1994). As to the GS equation, its 2D solu- 
tions with nonzero pressure were only obtained for pure 
hydrodynamical flows (see Beskin 1997 for review) when 
the isotropic pressure model is correct. 

On the other hand, for real astrophysical flows where 
the effects of finite pressure may play an important role, 
the isotropic pressure model does not seem to be sufficient. 
For example, in the solar wind where the free path ly ex- 
ceeds the Larmor radius rn by many order of magnitude, 
so that lu/rn ~ 10^ in the vicinity of the Earth, the effect 
of the anisotropic pressure may be important. 

Of course, it is necessary to remember that the ideal 
MHD approximation (including the double adiabatic ap- 
proach under consideration) is not good for real astrophys- 
ical objects, and so the dissipative processes can play an 
important role (Leer & Axford 1972; Hu, Esser & Habbal 
1997). Moreover, in the solar wind the stochastical com- 
ponent of the magnetic field is of the order of the regular 
one. Such fluctuations drastically change the velocity dis- 



tribution function. As a result, the pressure anisotropy 
P||/P^ Ki 1.2 is not actually large (see e.g. Marsch et al 
1982). May be, for this reason the effects of the anisotropic 
pressure on the outflow structure within the GS approach 
were analyzed in a small number of papers (Asseo & Beau- 
fils 1983; Tsikarishvih, Rogava & Tsikauri 1995), the GS 
equation itself being formulated only for the planar stat- 
ical configurations with zero longitudinal electric current 
(Notzel, Schindler & Birn 1985). 

In this paper we give the general equations describ- 
ing the axisymmetric stationary ideal MHD fiows with 
anisotropic pressure. Then, we shall present the results 
of our preliminary analysis of the effects of the anisotropic 
pressure (such an analysis has not been produced previ- 
ously either). 

2. BASIC EQUATIONS 

Let us consider a stationary axisymmetric (and nonrela- 
tivistic) MHD flow. Then the magnetic field can be written 
in the form 

B = — ^ e^. (1) 

ZTTtn wc 

Here '^{r,0) is the magnetic flux, m — rsinfl, and 
/(r, 9) is the total electric current flowing inside the region 
5" < ^'(r, 0). Assuming that the magnetosphere contains 
enough plasma to screen the longitudinal electric field E\\, 
one can write down 



E 



27rc 



(2) 



the 'field angular velocity' Vl-p = flp{'i>) being constant on 
the magnetic surfaces ^'(r, 0) = const. Finally, the conti- 
nuity equation V ■ (nu) = together with the freezing-in 
condition E + vxB/c = and the Maxwell equation 
V X E = result in 



V 
P 



B + rjpzue 



(3) 



where p is the mass density and 77 = '7(^) is the particle- 
to-magnetic flux ratio. Clearly, these expressions do not 
depend on the pressure anisotropy. 
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The next step is the transformation of the Euler equa- 
tion 



p(v • V)V = -V/3Pa/3 



^[VxB]xB + pg. (4) 



According to Chew, Goldberger & Low (1956), for colh- 
sionless plasma with an anisotropic pressure 



(5) 



it is possible to introduce two extra invariants on the mag- 
netic surfaces ^(j-, 0) — const (the so-called double adia- 
batic approximation resulting from the concervation of two 
adiabatic invariants) 



P B"^ 



El. 

pB- 



(6) 
(7) 



They correspond to polytropic equations of state with 
Fll = 3 and Fj^ = 1. As a result, integrating v and Lp com- 
ponents of equation (^), one can obtain the known expres- 
sions for the energy and the angular momentum (Asseo & 
Beaufils 1983; Tsikarishvili, Rogava & Tsikauri 1995): 



2 p 
L{^) = w^r sin6' - 

Here 



3 P, 
2 p 
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^-+2^(^-^^ 

(1 - /3a) . 

P - P 

^ n s 
P2 



GM 



, (8) 
(9) 



(10) 



is the anisotropic pressure parameter. Together with 
J7f(^) and ?/(^') these six invariants determine all the 
characteristics of a flow. In particular. 



= CTl 

27r 'l 



L — Qpw^ 



1 - M2 - /3a 



(11) 

(12) 



where = Anr]'^ / p is the Alfvenic Mach number. Thus, 
in the anisotropic case the Alfvenic singularity has the 
form 

A = 1 - - = 0. (13) 

It is necessary to stress that ( pi] ) gives an implicit ex- 
pression for the current / because / through B^ is con- 
tained in the rhs of this equation as well. Nevertheless, 
the standard procedure of the determination of flow pa- 
rameters in a given poloidal magnetic field Bp remains the 
same. The Bernoulli equation (g|), which can be rewritten 
in the form 



r(V*)2 



E-El-lEll 
p 2 p 



GM 
r 



(1 - NP - /3a)2 

-2n.^f^F(l-/3a)^^^^, (14) 



taken together with ([III ) and with the definitions (|^), (|^) 
determines implicitly the Mach number and /3a as 
functions of (or Bp) and all the six invariants 

^M\V^f,E,L,nF,Tl.si,S2\, (15) 

/3a = /3a[(V^')2, L, Qp, 77, Si, S2]. (16) 

It means that all the other parameters (/, v^, etc.) can 
be found from given values of Bp and of all the six in- 
variants as well. This fact is fully analogous to spherically 
symmetric flows (Bondi 1952; Parker 1958) where the flow 
structure is known and the algebraic expressions for the 
integrals of motion are enough to determine all the char- 
acteristics of a flow. 

On the other hand, if the flow is not spherically symmet- 
ric, the poloidal magnetic field is to be determined from 
the V^" component of equation (|^). After numerous but 
elementary calculations it can be presented in the follow- 
ing compact form 



1 



(1 - AP - (3^) V'=* 
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, 1 dsi , 1 dsp 



Here 



Si d* 



' S2 d^E* 



GAI 



n 



2 w^ nljl - /3a) - 2w^nFL{l - /3a) + AP 

1 - AP - 13, ■ 



(18) 



Vfc is a covariant derivative, and the operator d/d'^ acts 
on the invariants i?(*), IA^)-J^f{'^), and ri{^) only. Ac- 
cording to (||)-(0) and (p^-(|l^), this trans- field (stream, 
generalized GS) equation contains the unknown function 
^'(r, 0) and the six invariants only. As to the GS equation 
itself describing statical configurations, it can be obtained 



from (|17D in the limit fip ^ 0, 77 ^ (and hence A/I^ 
but rjL — > const. It gives for gravity- free case 



0), 



Vfc 



1 



(l-/3a) V'^* 



+ ^(l-/^a)/ 



d/ 
d* 



\ p 2 p 



1 dsi 



1 ds2 



^'Pn--^-^^^'Ps--d=^- (19) 



Si 



S2 d* 



Relations (H)-(|T^) remain true even for a more general 
bounded anisotropy model (Denton et al 1994) 



(vV) 
(vV) 



Pr,B^ 



P- 



El 

pB 



= 2Vt 



= -Vt 



P2 

1 

Tb 



(20) 
(21) 



where Vt is a nondissipative energy exchange term. In 
this case we have the concervation of the 'total entropy' 
S = S{^) 



(v-V)5 = i^(v.V)si 



pP(vV)s2=0, (22) 
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and so the generalized GS equation ( |l7| ) has a form 



Vfc 



.(l_M2-/3a)V'=* 



+ 



3271^ d fG_ 

-16^3 ^ Q ^23) 
d^f ^ ' 

Introducing the effective pressure Poff = ^cfr(^) as 

AP,^^p^(El + l^\.^S, (24) 

VP 2 p y 

which is equivalent to the ordinary thermodynamic rela- 
tion dP — pdw — pTds, one can obtain instead of (|9|) 



1 



^(l-/3a)V*--^' 



+167r3— ^ 
d^f 



d^- 
= 0. 



(25) 



For / = it has the form found by Notzel, Schindler & 
Birn (1985) for planar geometry. 

Equation ( p^ ) is a partial equation of the mixed type. 
Using the implicit relations resulting i n (n^ )"(p^), one can 
find for the second-order operator in (h 




- 4:7T- 



= 0, 

(26) 
(27) 



-3^ 



B^ 



in 



P ^2 

pv^ B^ 



3P — P \ 3P — P B^ 
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pi;2 



P2' 




(29) 



The form ( ^ ) of equation ( |l7| ) coincides with the one 
found for an isotropic pressure (see e.g. Sakurai 1990; 
Beskin 1997). Hence, one can conclude that the trans- 
field equation is elliptical for D > and D < —1, and 
hyperbolic for — 1 < Z) < 0. It means that the trans- field 
equation changes from elliptical to hyperbolic on the sur- 
faces D = (fast and slow magnetosonic singularities) and 
13 = —1 (cusp singularity). 

It is necessary to stress that, as in the isotropic case, 
D can be presented as £> — diA + d2B^/B^, so that 
D (X A for Bip = 0. This property reflects the fact that for 
k II B fast or slow magnetosonic velocity coincides with 



the Alfvenic one. In particular, for a nonrotating flow 
(r^F = 0, i = 0), when, according to ([Tl|), B^ — 0, we 
have 

n = (1 - m2 - /3a) (l - 3-^) . (30) 

Using now the definitions (0|) and (p7|), one can obtain for 
sonic, Alfvenic, and cusp velocities in anisotropic media 



3^, 
P 

P2 
A-Kp 

3— 

p Anp 



P ~ P 



P P 



p2 

— 

2 



B' ^Ps 

^ + 2 — 

Anp p 



(31) 
(32) 

(33) 



As we see, the velocities (|3l|)-(p3|) coincide exactly with 
the characteristic velocities propagating in anisotropic me- 
dia (Clemmow & Dougherty 1969). Recall that the con- 
dition v\ < corresponds to the firehose instability, and 
the condition i;2 < - to the mirror instability. 

3. TWO EXAMPLES 

3.1. Free outflow 

As a first example, let us consider a free (i.e. without 
gravity effects) outflow along a strong monopole magnetic 
field from the surface of a nonrotating sphere (with radius 
ro, magnetic field Bq, pressure Pq <C Bq/Att, mass density 
Po, and radial velocity vq). It is natural to assume that 
here Ps(ro) = Pn('ro) = Po- As a result, one can find for 
the radial velocity 



5 / 2r2 

2n^-5;^ 




where 



It gives 



Po 



-12p, 



div^vi) _ 



dr 



D 



where iV^ — AprQ/r'^, and 



D = 1 



Pn 

i — 

9 



(34) 



(35) 



(36) 



(37) 



which corresponds to the longitudinal sonic singularity in 
(p?!). On the other hand, the same as for isotropic media, 
for a nonrotating wind the magnetosonic surface coincides 
with the Alfvenic one, the singularity for 1 — M2 — /J^ = 
being absent. 

Comparing now expression (^^ for v{r) with the appro- 
priate one corresponding to an outflow with an isotropic 
pressure for the polytropic equation of state with F = 3 
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and the same energy E, one can conclude that the differ- 
ence only takes place for small enough pressure Pq <C PqVq 
when 



V, 



(iso) 



= 1 



V, 



(aniso) 



Vo 



2E 



2E 



(38) 
(39) 



For Pq :§> pqVq the limiting velocity does not depend on 
the pressure model: « 2E. Thus, the difference takes 
place for r ~ rg only. 

It is also interesting to determine how the anisotropic 
pressure affects the position of the Alfvenic surface rcr co- 
inciding for rip = with the fast magnetosonic one. As 



(40) 



it is easily check that j3^/M'^ <C 1 for both p <C 1 and 
p ^ 1, so that for Pq <C /OqUq and Pq ^ PqVq the position 
of the Alfvenic surface 



(41) 



does not depend on the pressure anisotropy. Finally, as 



Pn 
Ps 



,,2 ^2 



(42) 



the anisotropy of the pressure tensor increases with dis- 
tance from the central object. As has already been 
stressed, it can be realized only if dissipative processes 
(heat conductivity, pla sma instability, etc.) play no role. 
Relations (^) and ( p^ ) also mean that the conditions for 
the firehose and mirror instabilities are not satisfied up to 
the Alfvenic surface. On the other hand, the firehose insta- 
bility takes place in the supersonic region for r > rfii ^ rcr, 
where 



2 1 
r .— 



P \vo 



(43) 



3.2. Parker outflow 

As a next example we consider an outflow from the grav- 
ity center of mass M in a monopole magnetic field (Parker 
ejection). Analyzing the energy equation (H), we have 




(44) 



As a result, now iV„ = 4prQ/r^ — 2GM/w§r, and D is deter- 
mined by (37) again. Hence, for the position of the sonic 
surface we have 

''^=^''Gm'o- ^''^ 

On the other hand, one can easily check that equation 
di^ ) has a real solution for < ro only. This result is in 
agreement with a well-known fact that no transonic out- 
flow is possible for a gas with a polytropic index F = 3. 
Simultaneously, for Parker ejection the conditions of the 
firehose and mirror instabilities up to Alfvenic surface are 
not satisfied either. 

4. CONCLUSION 

Thus, we have generalized the GS equation for the flows 
with anisotropic pressure. It allows a quantitative analy- 
sis of the effect of anisotropic pressure upon the 2D out- 
flow structure within ideal double adiabatic approxima- 
tion. As has been shown on simple examples, the pres- 
sure anisotropy itself may be very large, but not in all 
cases does it result in a large disturbance of the flow pa- 
rameters. Glearly, the additional analysis (say, within the 
Weber-Davis model) is necessary to clarify all the details 
of the flow with anisotropic pressure. 
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